AN IDENTITY FOR SUMS OF POLYLOGARITHM FUNCTIONS 
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00 i Abstract. We derive an identity for certain linear combinations of polylogarithm functions with 

^^ I negative exponents, which impHes relations for linear combinations of Eulerian numbers. The coeffi- 

C^ i cients of our linear combinations are related to expanding moments of Satake parameters of holomor- 

^^ I phic cuspidal newforms in terms of the moments of the corresponding Fourier coefficients, which has 

S . applications in analyzing lower order terms in the behavior of zeros of L-functions near the central 

^ \ point. 
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1. Introduction 



Li,(a;) = ^k-'x^. (1.1) 
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\Q . If s is a negative integer, say s = —r, then the polylogarithm function converges for |x| < 1 and 

^^ i equals 

O , Y^r 

OO . , , 2^7=0 
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O . Li_.(a;) = ,^j;j,^, , (1.2) 

where the ( '. ) are the Eulerian numbers. The Eulerian number ( '^ ) is the number of permutations 



of {1, . . . , r} with j permutation ascents. One has 



"^ =i:(-irf't'V^'-^+i^ (1-3) 



e=o 
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From (|1.2I) we immediately deduce that, when s is a negative integer, Lis(x) is a rational function 
whose denominator is (1 — x)'"*'. Thus an appropriate integer linear combination of Lio(x) through 
Li_„(x) should be a simple rational function. In particular, we prove 

Theorem 1.1. Let a^^i be the coefficient ofk"^ in Y[i=o{^'^ ~ J^)' ^^d. let h^^i be the coefficient of k^ 

in {2k + 1) nj=o(^ ~ j)(^ + 1 + j)- Then for \x\ < 1 and d > 1 we have 

(2£)! x^(l + x) 






X 



i2£+l 



(2£+l)! 



X (1 + x) 

(l-x)2^+2' 



(1.5) 



We prove Theorem 1 1.1 1 in ^ While Theorem 1 1.1 1 only applies to linear combinations of polylog- 
arithm functions with s a negative integer, it is interesting to see how certain special combinations 
equal a very simple rational function. One application is to use this result to deduce relations among 
the Eulerian numbers (possibly by replacing x with 1 — x when expanding); another is of course 
to write Li_„(x) in terms of Li_„_|_i(x) through Lio(x). The coefficients a^^ and 6^^ which occur 
in our linear combinations also arise in expressions involving the Fourier coefficients of cuspidal 
newforms. We describe this connection in greater detail in ^ these expansions are related to un- 
derstanding the lower order terms in the behavior of zeros of L-functions of cuspidal newforms near 
the central point, (see nMiBII for a complete analysis). 



2. Proof of Theorem [UJ 
Before proving Theorem 1 1.1 1 we introduce some useful expressions. 
Definition 2.1. Let 



C2i 



l[{f-f) = (2£)!/2, C2,+i = {2i + l)l[{e-j){e + l+j) 

j=0 j=0 

Define constants Cm.r as follows: Cm,r = ifm ^ r mod 2, and 



{2i + l)\. (2.1) 
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(1) for r even, Co,o = 0, C2kfi = {—^Y'^for k >\, and for 1 < i < k set 






C2£ J-jJ C2£ {k-ey. 



(2) /or r o JJ and < i < k set 



C2£+l -^J^ C2£+l (fc-£)! 

A^o?e Cm,T- = ifm < r. Finally, set Br{x) = Ylm=o'^rn,ri—x)"^^'^ for \x\ < 1. Thus for r = 2i >2 
we have 



e-i 



i?2.(x) = E ^mM-^r^' = E ^^^^ n(^' - 3') 1 (-^)'- (2.4) 



, C2i 

m=0 fc=l \ j=0 

Immediately from the definition of Cr we have 

_ C2i _ C2i+1 

"''-' - T - 2£(2£ + 1) ' ^^-^^ 

as well as 

C2£+2 = (2£ + 2)(2£+l)c2,, C2,+3 = (2£ + 3) (2£ + 2)c2,+i- (2-6) 

While the definition of the Cm,r's above may seem arbitrary, these expressions arise in a very 
natural manner in number theory. See UMiBB for applications of these coefficients in understanding 
the behavior of zeros of GL(2) L-functions; we briefly discuss some of these relations in ^ 

Proof of Theorem [O] We first consider the case of r = 2£ even. We proceed by induction. We 
claim that 

^-(-) = ti^'^'f.f-^'))^-^^^ (2.7) 

for all L 

We consider the basis case, when (. = 1. Thus we must show for |x| < 1 that B2{x) = — s(l + 
x)/(l — x)^. As r = 2, the only non-zero terms are when m = 2k > is even. As C2 = 2 and 
C2fc,2 = i-l)^~^^k^ for A; > 1, we find that 

oo oo (-1 I \ 

B2{x) = Y.^-lf^'k\-xri' = -Y^k'x' = -U_2{x) = -^}^, (2.8) 

k=l k=l ^ ' 

which completes the proof of the basis step. For the inductive step, we assume 

and we must show the above holds with £ replaced by £ + 1. We apply the differential operator 

, N 2 

X—] -f (2.10) 

dx / 
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to both sides of (12.91) . After canceling the minus signs we obtain 

|:-.-,u('ri'(*'^-/^))-' = if^i^. <2->i) 

where the last line follows from (12.61) . which says C2^+2 = (2^ + 2)(2£ + l)c2^. Thus (12.71) is true 
for all i. 

As we have defined a£^i to be the coefficient of k^ in Y[j=oi^'^ ~ J^)^ (12.71) becomes 

j:j:a.,A:^x'= = C2,^^Ai±|l_. (2.12) 

A:=l J=0 ^ ^ 

The proof of Theorem 1 1.1 1 for r even is completed by noting that the left hand side above is just 

ae,2£^i-2£{x) H h a^,oLio(x). (2.13) 

The proof for r = 2£ + 1 odd proceeds similarly, the only significant difference is that now we 
apply the operator 

d \ f d 



X—] + a;— - £(£ + 1), (2.14) 

dx J \ dx J 

which will bring down a factor of [k — t}{k + 1 — tj. D 



3. Connections with number theory 

We now describe how our polylogarithm identity can be used to analyze zeros of L-functions 
near the central point. Katz and Samak UKaSa l conjecture that, in the limit as the conductors tend to 
infinity, the behavior of the normalized zeros near the central point agree with the N -^ oo scaling 
limit of the normalized eigenvalues near 1 of a subgroup of U{N) {N x N unitary matrices); see 
[IDMI EH |GUl EEl EMl ULSl iKaSal IMITTI iRol iRubl lYoJ for many examples. While the main terms 
for many families are the same as the conductors tend to infinity, a more careful analysis of the 
explicit formula allows us to isolate family dependent lower order terms. 

Our coefficients Cm,r are related to writing the moments of Satake parameters of certain GL(2) 
L-functions in terms of the moments of their Fourier coefficients, which we briefly review. Let 
Hl(N) be the set of all holomorphic cuspidal newforms of weight k and level N; see [.Iw2J for 
more details. Each / e Hl{N) has a Fourier expansion 

oo 

fi^) = ^o.fin)einz). (3.1) 

re=l 
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Let \f{n) = af{n)n~^''~^^^'^. These coefficients satisfy multiplicative relations, and |A/(p)| < 2. 
The L-function associated to / is 

71=1 P " / 

where xo is the principal character with modulus A^. We write 

\f{p) = af{p)+Pjip). (3.3) 

For p^ A^, af{p)(3f{p) = 1 and |a;/(p)| = 1. If p|A^ we take af{p) = Xf{p) and Pf{p) = 0. Letting 

denote the local factor at infinity, the completed L-function is 

A{s,f) = L^{s)L{s,f) = efAil-sJ), e/ = ±1. (3.5) 

The zeros of L-functions often encode arithmetic information, and their behavior is well-modeled 
by random matrix theory [CFKRSl IKaSai IKeSn3L The main tool in analyzing the behavior of 
these zeros is through an explicit formula, which relates sums of a test function at these zeros 
to sums of the Fourier transform of the test function at the primes, weighted by factors such as 
Q^/(p)™ + Pfip)"^- For example, if is an even Schwartz function, (f) its Fourier transform, and 
I + i'jf denotes a typical zero of A(s, /) for / G H^{N) (the Generalized Riemann Hypothesis 
asserts each 7/ G M), then the explicit formula is 

logA^' 



Yl SZ'^p'/ 



\Ht(N)\ ^^ ^^^ \'^ 2n 

_ -4(0) 1 ^ - af(^p)m ^ p^^p^n. i^gp g /^^ logP V .3 g. 

logiV \H;iN)\ ^^Z^^^ A. Z. pm/2 logAT '^ V logivj ' ^ ■ ^ 

see niLS[|Mil3ll for details and a definition of A(0) . Similar expansions hold for other families of L- 
functions. Information about the distribution of zeros in a family of L-functions (the left hand side 
above) is obtained by analyzing the prime sums weighted by the moments of the Satake parameters 
(on the right hand side). Thus it is important to be able to evaluate quantities such as 

t^^J^Mp^ + MpD (3.7) 

for various families of L-functions. 

For some problems it is convenient to rewrite a/(p)™ + /3/(p)™ in terms of a polynomial in 
Xf{p). This replaces moments of the Satake parameters a/(p) and Pf{p) with moments of the 
Fourier coefficients Xf{p), and for many problems the Fourier coefficients are more tractable; we 
give two examples. 

First, the p* coefficient of the L-function of the elliptic curve y"^ = x^ + Ax + _B is p^^/^ J2x mod p 
( ^ + ^+B ^ . j^gj-g ^^^ |g jj^g Legendre symbol, which is 1 if x is a non-zero square modulo p, if s = 
mod p, and —1 otherwise. Our sum equals the number of solutions to y'^ = x^ + Ax + B mod p, 
and thus these sums can be analyzed by using results on sums of Legendre symbols (see for example 
[lALMllMil2l ). 
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Second, the Petersson formula (see Corollary 2.10, Equation (2.58) of HILSII ') yields, for m,n > 1 
relatively prime to the level N, 

W^) 1^^ -.(/)A/(-)A/(n) =S^. + ((--)^/^^^^|§^) , (3.8) 

J k ^ ^ 

where 5mn = lifm = n and otherwise. Here the wji{f) are the harmonic weights 

wbW = Cn{2)/Z{1J) = C(2)/L(l,symV). (3.9) 

They are mildly varying, with (see HIwlllHLII ) 

N-'-^ <fe oonif) <fc iV-i+^; (3.10) 

if we allow ineffective constants we can replace A^^ with log A^ for A^ large. 

We can now see why our polylogarithm identity is useful. Using a f{p)+/3f{p) = Xp,af{p)f3f{p) = 
1 and \af{p)\ = \/3f{p)\ = 1, we find that 



afipf + Pfipf 


= MpY 


- 2 








c^fipr + Mp)' 


= ^fipf 


- 3A/(p) 








afipY + Pfipf 


= MpY 


- 4AKp)2 


+ 


2 




c^Apf + PApf 


= ^fipf 


- 5Xf{pr 
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5A/(p) 




c^f{pf + Pf{pf 


= ^fipY 


- Q^ApY 
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dXfipf - 


- 2 


c^ApY + PApY 


= ^ApY 


- 7>^ApY 
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UXfipf - 


- 7XAp) 


c^fipf + (3f{pr 


= ^ApY 


- 8Xf{pr 
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20Xf{pr - 


- IGXfipf + 2. 



(3.11) 



Writing af{p)™' + (3f{p)"^ as a polynomial in Xf{p), we find that 

m 
aApr + PAPr= E ^n^,rXApY, (3.12) 



r = m mod 2 



where the Cm,r are our coefficients from Definition [2T1 A key ingredient in the proof is noting that 

(1) C2k,2e. = C2k-i,2e-i - C2k-2,2i if ^ G {1, . . . , fc - 1} and /c > 2; 

(2) C2k+i,2e+i = C2fc,2£ — C2k-i,2e+i a i < k. 

We briefly describe the application of our identity, ignoring the book-keeping needed to deal with 



m <2. From the explicit formula (13.61) . we see we must understand sums such as 

where JF is a family of cuspidal newforms and Wr{J^) = Ylife^'^Rif) (^ simple Taylor series 
shows there is negligible contribution in replacing 0(mlogp/log_R) with 0(logj9/logi?)). As the 
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sums of powers of the Satake parameters are polynomials in A/(p), we may rewrite this as 

2^2^ 2^ pm/2 logR^^UogRj' '^ ^ 

p m=3 r^O ^ to \ to / 

r = 7r(. mod 2 

where ^^^^(p) is the r* moment of A/(p) in the family JF: 

A,,^(p) = :^^^ 5^ ^«(/)A/(p)^ (3.15) 

fes(p) 

We interchange the m and r sums (which is straightforward for p > 11, and follows by Abel 
summation for p < 7) and then apply our poly logarithm identity (Theorem ll.il) to rewrite the sum 
as 

^ ^ Ar,r{p)p''/^{p - 1) logp g / logp \ 

, ^ (p + ir+Mogi? ^{logRj- 

For many families we either know or conjecture a distribution for the (weighted) Fourier coef- 
ficients. If this were the case, then we could replace the Ar^jr{p) with the r* moment. In many 
applications (for example, using the Petersson formula for families of cuspidal newforms of fixed 
weight and square-free level tending to infinity) we know the moments up to a negligible correction 
(the distribution is often known or conjectured to be Sato-Tate, unless we are looking at families 
of elliptic curves with complex multiplication, where the distribution is known and slightly more 
complicated). Simple algebra yields 

Lemma 3.1. Assume for r > 3 that 

ArAP) = { ^'"^^ (3- IV) 

[ otherwise, 

and that there is a nice function qm such that 

oo 

gM{x) = M2X^ + Msx^ + ■ ■ ■ = J^MiX^- (3.18) 

e=2 
Then the contribution from the r > 3 terms in the explicit formula is 

Thus we can use our polylogarithm identity to rewrite the sums arising in the explicit formula 
in a very compact way which emphasizes properties of the known or conjectured distribution of 
the Fourier coefficients. One application of this is in analyzing the behavior of the zeros of L- 
functions near the central point. Many investigations have shown that, for numerous families, as 
the conductors tend to infinity the behavior of these zeros is the same as the N -^ oo scaling limit 
of eigenvalues near 1 of subgroups of U{N). 

Most of these studies only examine the main term, showing agreement in the limit with random 
matrix theory (the scaling limits of eigenvalues of U{N)). In particular, all one-parameter families 
of elliptic curves over Q(T) with the same rank and same limiting distribution of signs of functional 
equation have the same main term for the behavior of their zeros. What is unsatisfying about this 
is that the arithmetic of the families is not seen; this is remedied, however, by studying the lower 
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order terms in the 1-level density. There we do break the universality and see arithmetic dependent 
terms. In particular, our formula shows that we have different answers for families of elliptic curves 
with and without complex multiplication (as these two cases have different densities for the Fourier 
coefficients). 

These lower order differences, which reflect the arithmetic structure of the family, are quite im- 
portant. While the behavior of many properties of zeros of L-functions of height T are well-modeled 
by the A^ — > oo scaling limits of eigenvalues of a classical compact group, better agreement (tak- 
ing into account lower order terms) is given by studying matrices of size N = (logT)/27r (see 
HKeSnll |KeSn2l |KeSn3ll ). Recently it has been observed that even better agreement is obtained 
by replacing N with N^s, where A^cff is chosen so that the main and first lower order terms match 
(see [BBLM[ IDHKMSII ). Thus one consequence of our work is in deriving a tractable formula to 
identify the lower order correction terms, which results in an improved model for the behavior of 
the zeros. 
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